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ABSTRACT 

The SC and CS stars are thermal-pulsing AGB stars with C/O ratio close to unity. Within 
this small group, the Mira variable BH Cm recently evolved from spectral type SC (showing 
ZrO bands) to CS (showing weak C2). Wavelet analysis shows that the spectral evolution was 
accompanied by a dramatic period increase, from 420 to 540 days, indicating an expanding 
radius. The pulsation amplitude also increased. Old photographic plates are used to establish 
that the period before 1940 was around 490 days. Chemical models indicate that the spectral 
changes were caused by a decrease in stellar temperature, related to the increasing radius. 
There is no evidence for a change in C/O ratio. The evolution in BH Cm is unlikely to be 
related to an on-going thermal pulse. Periods of the other SC and CS stars, including nine new 
periods, are determined. A second SC star, LX Cyg, also shows evidence for a large increase 
in period, and one further star shows a period inconsistent with a previous determination. Mira 
periods may be intrinsically unstable for C/O « 1; possibly because of a feedback between 
the molecular opacities, pulsation amplitude, and period. LRS spectra of 6 SC stars suggest 
a feature at A > 15/im, which resembles one recently attributed to the iron-sulfide troilite. 
Chemical models predict a large abundance of FeS in SC stars, in agreement with the proposed 
association. 

Key words: stars: individual: BH Cm - stars: AGB and post-AGB - stars: spectral classifi- 
cation 



1 INTRODUCTION 

The observed properties of stars on the Asymptotic Giant Branch 
(AGB) are largely determined by the ratio of carbon to oxygen in 
their atmospheres. Carbon abundances in AGB stars increase due to 
dredge-up following their regular thermal pulses. After a series of 
these events, which happen every 10 4 -10 5 yr, the enhanced carbon 
may cause a transition from an oxygen-ric h star to a ca rbon-rich 
star, via the intermediate S-stars (e.g., iLlovd Evanil984l) . 

The SC stars iKeenan & Boeshaa3rFJ8(3) form a continu- 
ous spectral sequ ence intermediate between the S and C stars 
iCatchoole & Feast 19 7 H). They show very strong sodium D-lines, 
and strong CN bands Tstephensonlll973t) . and either weak ZrO 
bands (SC stars) or C2 bands (CS stars). The molecular abundances 
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indicate a C/O number ratio very close to unity, so that CO forma- 
tion leaves little C or O for the other molecules. 

The SC star BH Cruris has simultaneously been found 
to show a unique combination o f a lengthening period 
jBateson. Mcintosh & Venimord Il988t I Walker & Marine! 1 199 it 
Walker, Ives & Williams 1995), and evolution in spectral type from 
SC to CS iLlovd Evansl 19851) . with ZrO bands iCatchpole & Feast! 
119711) disappearing and C2 band appearing. The spectral evolution 
has b een interpreted as due to an increase in C/O ratio IWhitelockl 
1999). Both the period and abundance changes have been sug- 
geste d to be caused by a recent thermal pulse iWhitelock 1999; 
Wood & Zarro 1981). BH Cm could therefore present a unique case 
of real-time AGB evolution. 

In this paper we present an analysis of the light curve of 
BH Cm, to quantify the period and amplitude variation. Chemi- 
cal modelling is used to investigate whether the change in relative 
molecular abundances can be explained as due to a decrease in ef- 
fective temperature, at constant C/O ratio. We also present new pe- 
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riods for other SC stars, including two stars which may show period 
changes similar to BH Cm. 

2 PULSATION IN BH CRU 
2.1 The star 

BH C m (also known as He 120) was reported to be variable by 
Welch jl969|), and assigned a sp ectral type of SC 4.5-7 iKeenan 
1971; Keenan & Boeshaar 1980). The star is bright, reaching a 
visual magnitude of 7, with an amplitude between l m and 3 m . 
IWhitelock etal] l2000l) suggest a distance of d = 900 pc, from the 
Mira PL relation. At Galactic m etallicity, carbon sta rs form for ini- 
tial masses ^ 2 Mq iGroeneweeen & de Jond 19931 their Table 5). 
The galactic z-height of BH Cru of about 100 pc is consistent with 
such a progenitor star. A strong Sr 1 4607 A line shows evidence for 
s-process element enhancements. 



Optical photometry fr omlWalkedl 19791) and infrared photom- 
etry from White lock et alJ |2000) are listed in TableQ The colours 
are typical for unobscured carbon stars, and the IRAS flux densi- 
ties are consistent with photospheric emission. The K— [12]= 0.96 
colour gives a black-body temperature of 3000 K, identical to the 
stellar temperature ( section l3~3l below) . There is no evidence for on- 
going mass loss, neither from self-obscuration in the near-infrared, 
nor from a mid-infrared dust excess. The lack of mass loss is un- 
usual for such a long-period Mira. 



2.2 Recent period evolution 

We have analysed visual observations collected by the Royal Astro- 
nomical Society of New Zealand (RASNZ). Only data from indi- 
vidual observers contributing 15-20 observations or more for each 
star were used. The light curve was binned into 20-day averages 
and analysed using wavelet transforms, which have the advantage 
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Table 1. Visual and infrared m agnitudes of BH Cru, as given in IWalkeil 
Il979tl and lWhitelock et alJ<200of . 
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of being sensitive to chang es of the pulsation properties (period, 
amplitude) over time (e.g., Bed ding et alj l998>. We have used the 
weighted wavelet Z- transform (WWZ) (Foster 1996) developed at 
the AAVSO specifically for unevenly sampled data. 

The visual light curve is shown in the top panel of Fig. Q 
During the 1970s, the maxima tended to be flat and long-lived, 
with pronounced but shorter-lasting minima. From 1990 onward, 
the maxima became more pronounced and slightly brighter, with a 
hump on the ascending part that is very common among Mira vari- 
ables. In the first few years, the light curve also showed a secondary 
minimum superposed on the plateau of maximum. Double max- 
ima and minima of alternating depth are not uncommon amo ng the 
longer-pe riod Miras, e.g. R Cen (Hawkins, Matte i & FosteJ200lh . 
In R Dor l B edding etalJl998t) . a split maximum preceded a change 
in pulsation mode. 

The wavelet plot of BH Cru confirms the period evolution first 
reported by Bate son et alj, ll9881. The second panel of Fig.Qshows 
the pulsation frequency as function of time, where the contours rep- 
resent the significance level. The 'best' period is plotted in the bot- 
tom panel. From 1970 until around 1995, we see a steady increase 
in the period. Since 1995, the period has been more or less con- 
stant. The period was 425 days at the time of discovery (1969), 
500 days in 1990 and 535 days in 1999. This large increase of over 
25% within (at most) 20 years is unparallelled: among known Mi- 
ras, only T UMi <Gal & Szatmar\ll995l:lMattei & FosteJ 19951) has 
show n a comparable, but o pposite, change. The new case of LX 
Cyg iTempleton et all2 003) is discussed further below. 

Although the period evolution shows some possible fluctua- 
tions, it can be well fitted with a constant rate of change of 5 days 
per year, corresponding to 1.4% per cycle, between 1975 and 1995. 
(The pre-1970 part of the curve shows a reverse change, but this 
may be affected by edge effects in the wavelet analysis.) 

The semi-amplitude (third panel of Fig0 shows similar strong 
evolution. It was around 0.8 mag before 1980 and, in parallel with 
the period increase, the semi-amplitude rapidly increased to 1.4 
mag in about 1992. Except for a single faint maximum in AD 2000, 
the amplitude has since remained constant. 

Fig.|2|illustrates the correlation between the increases in am- 
plitude and period. It shows a roughly linear correlation, with 
the semi-amplitude increasing by about 50 mmag/cycle, or about 
5% per cycle. The correlation between period and semi-amp litude 
evolu tion mirrors that seen in R Hva IZiil stra. Be dding & Matteil 
120021). and also i n R Aql and S Ori jBedding et all2000h and Y Per 
iKissetalli oOO). 



2.3 Older period data 

Although BH Cru was only discovered to be variable in 1969, the 
region has been extensively covered by older surveys. But we have 
not found published records of BH Cru prior to 1969. The star is 
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Figure 2. The visual semi-amplitude, in magnitudes, versus the period in 
days, since JD 2441000 (AD 1970). The dotted line shows the lower limit 
of the Mira classification. 



absent from the CP catalogue (Gill & Kaptevn 1900) but this has a 
completeness limit of 9.2 mag and the blue magnitude of BH Cru 
is always below this. The deeper CD catalogue lThomell932T) goes 
to 10 mag, but also does not include it. These limits do not allow 
for an unambiguous interpretation. We have therefore searched the 
Harvard plates for evidence of its pre-discovery behaviour. 

The Harvard Plate Stack contain 96 good observations of BH 
Cru, 15 detections with uncertain magnitude, and 37 plates where 
the star is not visible. The plates are almost exclusively blue. Com- 
bined with BH Cru's B — V, the star exists near the limit of the 
plates while at maximum. The data are shown in Fig. [3] We per- 
formed a wavelet analysis using the 96 good and 15 uncertain deter- 
minations, and also included upper limits where they were fainter 
than 1 1.8mag. (The brighter limits do not constrain the variability.) 

BH Cru was a pulsating variable well before its discovery as 
such. Strong pulsations are evident, with an amplitude between 1 
and 2 magnitude. The plates before 1910 are consistent with little or 
no fluctuations, but the sampling is poor and this is not conclusive. 
The data are too patchy for a well-determined period, however a 
period of about 490 days is indicated around 1940. The period at 
1930 may have been slightly shorter (about 440 days). 

Compared to Fig.Q the period around 1940 was about 10 per 
cent longer than in the early 1970's. The data suggest that the period 
was not constant even in this early phase and there may have been a 
further increase shortly after 1940. Further data would be desirable, 
but the conclusion from the available data is that between 1945 and 
1970 the period must have decreased. The recent period instability 
may be part of a recurrent behaviour. 
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Figure 3. Photographic observations of BH Cm in 20-day bins. Upper 
limits, where a comparison star of the that magnitude could be seen, are in- 
cluded if at m=11.8, but brighter limits are not. The systematic photometric 
accuracy of the data is no better than 0.5 mag. Period and amplitude before 
JD2425000 have very low confidence. 



2.4 Physical changes in BH Cru 

The pulsation equation for Mira variables is 
log P = 1.5 log R — 0.5 log M + log Q, 



for first overtone pulsators, where the pulsation constant Q 
iFox & WooJl982l) . or 



(1) 
0.04 



log P = 1.949 log R — 0.9 log M — 2.07 (2) 

for fundamental mode pulsators lWoodll990h . P is the period in 
days and R and M are the radius and mass in solar units. 

In Mira variables, M decreases on times of scales of 10 - 
10 7 yr because of the stellar wind, but over a century it is approxi- 
mately constant. In any case, BH Cru shows no evidence for large 
mass loss. Q can be assumed as constant on this short time scale. 
The radius R is not directly known; it is however the only parame- 
ter whose changes can explain the period evolution of BH Cru. 

Miras are thought to be fundamental-mode pulsators, although 
this is not as well determined for the very long period variables. 
However, the observed optical radii for Miras are much larger than 
those predicted from the fundamental mode, and agree better with 
first overtone pulsators 1 Haniff, Scholz & Tuthill ll995l) . The dis- 
crepa ncy may be due t o missing opacity in the atmosphere models. 
Ivan Belle et alj {2(302) find increased diameters also at K, which 
they attribute to circumstellar (water) emission. Here we prefer Eq. 
Q since it better fits the atmospheres. 

The period of BH Cru increased by 25% since 1970, cor- 
responding to a 16% change in the radius from Eq. 0. For a 
stable bolometric magnitude, this would decrease the effective 
temperatur e by 7.5%: ~ 250 K at the temperature of 3200 K 
iLoidl. Lancon & Jeirgensenl200ll) . Fundamental-mode models re- 



quire a smaller radius change (12%) and temperature change (~ 
200 K). 

If the luminosity also increased, a smaller temperature de- 
crease would be required. The mean magnitude of the infrared light 
curve changed by no more than 0.1 mag, between 1980 and 1995, 
and the J— K colour reddened by no more than 0.05 mag. The bolo- 
metric correction at K therefore decreased by less than 0.05 mag, 
using Fig. 16 of IWhitelock et alJ feOOOh . The mean magnitude of 
the optical light curve was constant to within 0.5 mag but it is 
much more sensitive to temperature variations and varying molecu- 
lar opacity. Overall, the data is consistent with evolution at constant 
luminosity; any increase in the luminosity did not exceed 10%. 

Assuming an initial radius of 300 Rq, BH Cru will have ex- 
panded by 6 x 10 10 m, at an average velocity of ~ 100 m s _1 . Typ- 
ical small-amplitude Mira variables expand at 3-10 kms -1 during 
the pulsation cycle. The small average expansion velocity would 
therefore not be expected to significantly disrupt the pulsational 
behaviour. 



3 EQUILIBRIUM CHEMISTRY: ZrO AND C 2 
3.1 Spectral changes 

The spectral evolution of BH Cru is remarkable. Before 
1973, BH Cru showed e vidence for ZrO, but lacked C2 lines 
ICatchpole & Feast! Il97lh: the ZrO line was weak or absent 
around 1970 lKeenarll97ll) but strong in 1967 fctephensorll973l) : 
ICatchpole & FeastHl97lf reported it as weakly present on an un- 
dated plate. In 1980, line s of C2 had appeare d and the ZrO lines 
were no longer present jLlovd Evansl Il985l) . The spectra taken 
around 1980 cover a substantial part of the light c urve: cy clical 
variab ility cannot account for the observed evolution iLlovd Evansl 
1985). The few spectra taken around 1970 show ZrO lines of vari- 
able strength: during some part of the pulsation cycle, the pre-1973 
spectra may have been lacking in ZrO as well. But the strong ZrO 
bands visible at times before 1973 did not re-occur over any part of 
te pulsation cycle in 1980. 

The change has been interpreted as tracing an increase in the 
C/O ratio from < 1 to > 1 1 Whitelockl ll999h . associated with the 
dredge-up of carbon produced in a helium shell flash (a thermal 
pulse). However, below we investigate whether the change in effec- 
tive temperature of the star can also explain the abundance changes, 
without requiring a change in the C/O ratio. We note that a possible 
recurrent shift in spectrum between the presence of ZrO and C2 has 
been seen in TT Cen 1 Stephenson ! 1 9731) . 

A temperature effect is indicated by the fact that SC stars tend 
to have earlier subtypes than CS stars. Only two SC stars (see Ap- 
pendix^ are known to reach types later than SC5 (BH Cru and LX 
Cyg), while three of the five CS stars are classed as C7 or later. The 
temperature zero point for C and S stars subclasses differ: C5 cor- 
responds in temperature to about M0. However, SC7 and C7 have 
similar temperatures (Keenan & Boeshaar 1980), and the late CS 
stars have therefore lower stellar temperatures than any of the SC 
stars with known subclass. 



3.2 Temperature dependence 

To investigate the temperature and C/O dependence in isolation 
of other parameters, we model the photospheric composition as a 
molecular equilibrium (e.g.. lTsuiilll973l: IWvckoff & Clegdll978t 
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Figure 4. The equilibrium partial pressure of C2 (log), as function of tem- 
perature, for seven different values of the C/O ratio: 1.01, 1.05, 1.1, 1.2, 1.3, 
1.4 and 1.5 (lower curve to upper curve respectively) 
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Figure 5. Same as Fig.|4] but for ZrO 



Figure 6. Same as Fig.|4] but for YO 
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Figure 7. Same as Fig.|4] but for H2. Molecular hydrogen is so dominant 
in the model that the variation with C/O ratio is negligible. 



Tarafd ailll987t Isharo & Huebnetll 199(1) . We calculate the molecu- 
lar equilibrium composition by a direct minimisation of the Gibbs 
free energy of the system usin g a steepest descent technique 
( Whit e. Johnson & Dantzi3ll958l) . The free energy of the system 
is given by 

G = ^f iXi , (3) 

i 

where Xi is the number of moles of molecule i and fi is the chem- 
ical potential of molecule i, 

^(^) l +lnP+ln (D- (4) 

Here, (G/RT) i is the Gibbs free energy of formation of 
molecule i, P is the total pressure of the system (a model pa- 
rameter), and x — Xi. The Gibbs free energy of formation of 
a molecule is a temperature-dependent thermo-chemical property 
that c an be found for example in the JANAF tables 1 Cha se et alJ 
Il985l) . To estimate the value at other tempera tures, we used an in- 
terpolation formula. SharD & Huebner 1 1990) provide coefficients 
to a polynomial fit of the form 

AG = aT~ 1 + b + cT + dT 2 + eT 3 , (5) 



to the Gibbs free energies for 210 species using JANAF data. We 
added further species appropriate to AGB stars i n the same way. 

At chemical equilibrium, G is a minimum rZema nsk^ll957l) . 
The equilibrium composition is therefore the set of mole num- 
bers X — {xi} i ^ 1 which minimises G. Full details of the min- 
imisation technique subject to these constraints can be found in 
lMarkwickli200bT) . 

We adopt solar elemental abundances JAnders & Grevessel 
1989; Cameron 1973) except for carbon which is enhanced to vary 
the C/O ratio. Models were run for values of the C/O ratio, be- 
tween 1.01 and 1.5. The resulting equilibrium abundances of the 
molecules C2, ZrO and YO for a variety of temperatures and C/O 
ratios are shown in Figures [4j|6] A plot of the H2 partial pressure 
is given in Fig. 0to allow conversion of the partial pressures into 
fractional abundances. 

Figure [4] shows clearly that the photospheric abundance of 
C2 increases with decreasing temperature over the range 4000- 
2800 K. Furthermore, the model shows that oxides such as ZrO 
can still exist at C/O ratio > 1, although at reduced abundance 
(Fig. The relative abundances shown in the figure do not take 
the s-p rocess enhanc ement s into account. In BH Cru, [Zr/H] = 1.2 
lAbia & Wallersteinl 1 1998) : assuming a linear scaling, the ZrO 
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and C2 relative abundances are almost equal for C/0= 1.01 at 
T = 3300 K. lAbia & Wallersteirl < 19981) give C/O = 1.02. 

The photospheric temperature of BH Cru is taken as 3200 K 
(Loidl et al. 2001; these authors found a temperature close to min- 
imum of about 2800 K). At these temperatures, Fig. |4| shows that 
even a relatively small decrease in temperature leads to a large in- 
crease in the abundance of C2: a decrease of 300 K can yield a 
ten-fold increase. The peak abundance is reached at T ~ 2800 K. 

The disappearance of ZrO and appearance of C2 took place 
between 1973 and 1980 lLlovd Evans] 1 19851) . Assuming a constant 
rate of period evolution, the period increased by ~ 50 days or 
roughly 12% over this time. This corresponds to a radius increase of 
8% and a temperature decrease of 4% if the luminosity remained 
constant (assuming first overtone pulsation). The AT = 130 K 
for a current temperature of 3200 K would have increased the C2 
abundance by a factor of 3-5 and decreased the ZrO abundance 
by a similar factor. The YO abundance should have decreased as 
well (see Fig.EH, but this m olecule has only been weakly detected 
iKeenan &"B oeshaar 1980) and there is no observational mention 
of a change with time. 

3.3 Atmosphere models 

The model above has two shortcomings; we do not consider the ex- 
tended atmosphere, and no time-dependent chemistry is included. 
Chemical equilibrium is unlikely to hold at the conditions and dy- 
namic time scales within an extended pulsation-driven atmosphere. 
However, we do not attempt to address time-dependent chemistry 
here. 

The extended atm osphere is describ ed using the hydrostatic 
atmosphere models of Loidl et alj 120011) . These models include 
equilibrium chemistry for a large number of molecules (but exclud- 
ing e.g., Mg and Fe molecules and CaCl). Opacities are included 
for 7 carbon-star molecules only: CO, CN, CH, C 2 , HCN, C 2 H 2 
and C3. C/0= 1.01 is used. The models are shown in Fig.|8|for 
two different temperatures, corresponding to two different phases 
of the pulsation cycle. Comparison between the two temperatures 
confirms the strong dependence of the oxides and C2 on tempera- 
ture. At the higher temperatures, C2 is mainly present near the pho- 
tosphere. At the lower temperature, the abundance in the extended 
atmosphere is much higher. Model spectra at the two temperatures 
are shown in Fig. [9] 

3.4 Optical spectra 

The spectra of oxygen-rich Miras variables show few atomic lines 
due to the absorbing molecular blanket which hides the photo- 
sphere. In SC stars, this blanket is transparent and the spectrum 
is covered in lines tTsuii 1 964T) . 

Infared spectra of BH Cru are presented bv lLoidletallfeOOll) : 
they show strong CN and weak C2 bands. We obtained optical spec- 
tra of BH Cru and some comparison SC/CS stars using the double- 
beam spectrograph at the MSSSO 2.3-m telescope, centred at the 
6500A ZrO band. The data were taken on 2003 May 18 and 21, 
under non-photometric conditions. The dispersion is 0.55A/pixel 
and the resolution is about lA. The spectra are flat-fielded using an 
internal lamp and sky-subtracted. We did not correct for the atmo- 
spheric and instrumental response. 

The spectra are shown in Fig.^| Typical noise per wavelength 
pixel is 1%. The many lines repeat b etween the spectra, and ar e 
largely due to s-enhanced metals (e.g. lAbia & W allerstein 19981) . 



The 6708A line is present in all other spectra but weak or absent in 
BH Cru. In all cases it is weaker than the adjacent Lai 6710A. The 
lithium line o verlaps w ith another line and the identification is not 
unambiguous iRevniers et al. 2002). 

The model spectrum of BH Cru (see previous section) is 
shown in the top diagram. The model is for 3200 K: the 2900 K 
model also gives an acceptable fit. The good fit obtained around 
6500A confirms the absence of the ZrO band (6474-6540A) in BH 
Cru: the good fit is obtained in spite of this band not being included 
in the opacities in the model. All other stars which were observed 
show the band, although with varying strength. The broad 6350A 
ZrO band is similarly seen in all stars except BH Cru. 

BH Cru shows possible absorption lines at 6053 & 6070A, 
which differ from the other stars and may be from the C2 Swan- 
band. Other lines of the Swan bands at 6190 & 6480A are not 
detected, but the confusion at our resolution is too high to rule 
out their presence . C2 is present in BH Cru longward of 7000A 
iLoidl et allioOlh . 

The double CaCl band at 6200A is seen in TT Cen only. This 
band is common in carbon stars, but is often blended and masked 
by strong CN bands. The CaCl abundance depends very little on 
the C/O ratio (Fig.|l 1>. The CaCl band strength tends to vary over 
the pulsation cycle JSanford 1950), being strongest near minimum 
light: TT Cen was observed near minimum. TT Cen may be com- 
pared to VX Aql iWvckoff & WehingeJl976T) . in which CaCl also 
is the strongest molecular band . Both these stars are long-period 
Miras. Clegg & Wvckoff ( 1977) show that CaCl is favoured by high 
density environments. 

3.5 Dust 

Stars with C/O close to unity have difficulty forming dust: the re- 
moval of C and O from the chemistry means that neither silicates 
nor carbon grains can form. It is therefore not surprising that the 
SC stars show little evidence for significant dust emission. A few 
ISO spectra have very poor S/N at mid-infrared wavelengths, and 
weak features would not be seen. 

We therefore retrieved the IRAS LRS spectra of the seven 
SC/CS stars with \2-ym fluxes above 15 Jy: BH Cru, AM Cen, FU 
Mon, UY Cen, S Lyr, VX Aql, and WZ Cas. The spectra mostly 
show adequate but low S/N (the brightest star in the sample, UY 
Cen, only has F\2 = 55 Jy). In order to emphasize any excess, we 
have subtracted a pseudo-continuum represented by a A -2 power 
law that was scaled to the flux at 14 /im for each star. 

The top panel of Fig. ll2l shows the resulting spectrum of WZ 
Cas. The impression of a weak silicate feature is deceptive: it re- 
flects a deep CS absorption feature at 8^im. 

None of the other SC LRS stars resembles WZ Cas. We have 
co-added the remaining six spectra and again subtracted a pseudo- 
continuum to produce a 'typical' SC-star LRS spectrum (lower 
panel of Fig. fL?! . Compared to WZ Cas, the continuum-subtracted 
spectrum rises beyond 15 /im. The wavelength coverage is insuffi- 
cient to decide whether this is a dust continuum or a broad feature. 
However, the rise resembles that seen in some young stellar objects, 
where it is part of a feature peaking at 23/^m. This fe ature has pre- 
viously been attributed to FeO, but Keller et al. 1 2002) found a fit to 
the iron-sulfide troilite. lHonv et alj 120 02) report the same feature 
and presumably the same mineral in planetary nebulae, and sug- 
gest a carbon-star origin. If the SC stars show the same feature, it 
would strongly support the identification with troilite, as formation 
of this mineral should benefit from the removal of O from the dust 
formation. 
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Figure 8. Results of the atmosphere models. Top-left: temperature, density and opacity for a stellar temperature of 3200 K. The stellar surface is on the right. 
The opacity scale runs from to 1. Top-right: same, for a stellar temperature of 2900 K. Bottom: Molecular partial pressures for the two temperatures. 



Figures fT3l and l 14l show the calculated equilibrium abundances 
of FeO and FeS. As expected, FeO declines sharply as the C/O ratio 
increases. FeS also shows a decline but is less sensitive especially 
at photospheric temperatures. At all but the lowest temperatures, 
FeS is the more abundant molecule by three orders of magnitude. 
The formation of an iron sulfide such as troilite would therefore not 
be unexpected. However, the evidence for a troilite feature is weak 
and this point re mains to be confirmed. 

iFerrarotti et alJ J2000) suggest that FeSi would be the main 
dust component in S stars. This mineral has bands at 32 and 50 
micron. ISO SWS spectra are available for three SC stars: WZ Cas, 
S Lyr and W Cas. The S/N is low and no emission can be seen 
longward of 30 micron. 

Evidence for mass loss in BH Cru comes from a blue-shifted 



(by 19kms _1 ) possible rubidium line lAbia & Wallersteinl il99^) . 
Otherwise, the SC stars show little evidence for significant mass 
loss. (Of the stars in Table I All only FU Mon was detected in CO 
and a mass-loss rate of 10 -7 yr" 1 derived with a v ery low ex- 
pansion velocity (2.8 kms -1 ) Jjorissen & K napp 1998).) This is in 
marked contrast to the oxygen-rich Miras and carbon Miras, even 
though the SC stars form an intermediate evolution. It is however 
understandable, as the mass loss is driven by radiation pressure on 
dust. S C stars are expected to form a phase of mass-loss interrup- 
tion JZiilstra et alll992l) while evolving from an oxygen-rich to a 
carbon-rich star. 
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Figure 9. The model spectrum of BH Cm for two temperatures: 3200 K (drawn) and 2900 K (dashed) 
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Figure 11. Same as Fig.E| but for CaCI 



4 LX CYG: ANOTHER BH CRU 

Large period changes similar to BH Cru are very rare among Mira 
variables: IZiilstra & Beddiri3 i200 3) describe 6 cases, including R 
Hya. They estimate that the type of evolution seen in BH Cru oc- 
curs in ~ 1 per cent of well observed Miras, over a century of 
observations. To find one on-going case among the small group of 
SC/CS stars (16 in total) would be unexpected. To find a second 
one would suggest a correlation between the rare spectral class and 
the period evolution. 

There is in fact a second case of period evolution among the 
SC stars. LX Cyg is an SC3 star with a catalogued (GCVS) period 
of 465 days. An increase of its p eriod was independently discov- 
ered bvlAlksne & Alksnisl 1 19851) and by iBroens et alj fcOOd) and 
lTempletoti^^lT]2003TriTempletori et aljl2003l) have analyzed an 
extended data set. They find that historical records of LX Cyg show 
the period was fairly stable at about 460d from its first measurement 
around JD 2426000 to 2439000, after which it increased to around 
580 days. We note that BH Cru and LX Cyg are the two SC stars in 
Table lAll with the latest subtypes and longest periods. A spectrum 
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Figure 12. IRAS LRS spectra of SC/CS stars, from which a A -2 pseudo- 
continuum has been subtracted. Top panel: WZ Cas, the only star showing 
a silicate-like feature. Bottom panel: average of six other SC stars. 



kindly taken for us by James Bryan in 2003 confirms that LX Cyg is 
an SC star, but does not allow a current subclass to be determined. 

We carried out a wavelet analysis of the AAVSO data, identi- 
cal to the one shown for BH Cru. The result is shown in Fig. \^\ 
the axis scales are the same as in Fig.Q to facilitate the compari- 
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Figure 10. Optical spectra of seven SC/CS stars. The inset shows the (saturated) Ha emission line of BH Cm. The spectra are not response-corrected or 
flux-calibrated. The strong band at 6857A is telluric. The spectrum of BH Cm is shown together with a model for 3200 K. The y-axis gives total counts 
X 10~ 4 per wavelength pixel; the noise is negligible, and all line features are real. A few lines/bands are labeled; indicated but unlabeled lines/blends are due 
toCN. 
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Figure 14. Same as Fig.|4] but for FeS 

son. The rate of period change in LX Cyg is very similar to that in 
BH Cm. The pulsation amplitude is somewhat larger and has not 
changed significantly, and the total change in period is larger. 

A redetermination, or in several cases a first determination, 
of the periods of the SC/CS stars is given Appendix There we 
show that our period of VX Aql (609 days; unknown subtype) is 
inconsistent with a earlier determination (470 days). We have in- 
sufficient data to confirm a period cange, but it seems possible that 
VX Aql is a third case of a SC star with an evolving period. If true, 
all SC star with periods longer than about 470 days show this ef- 
fect. Many n ormal Miras with such periods also show fluctuations 
(Ziilstra & Bedding 2003), but limited to about 10 per cent in pe- 
riod. The SC stars show much larger changes. Even without VX 
Aql, a link between the C/O ratio near unity, and the unstable peri- 
ods, seems likely. 



5 DISCUSSION 

A few Miras are known to show changing peri ods. T he well-known 
case of R Hya was studied by IZiilstra et alj 120021) who find evi- 
dence for a stable period of 495 days before AD 1800, followed by 
a declining period, and a newly stable period of 385 dats since AD 
1950. R Aql shows similar behaviour. TY Cas was recently added 



to the list of variable periods lHazen & Mattell2003l). But overall, 
period evolution among Miras is rare. iGreaves^^Howarfhl l2000l) 
examined a sample of 100 Mir as, and apart from the prev iously 
known case of T UMi ( Gal & Szatmarv 1995; Mattei & Foster 
1995; Szatmarv et al. 2003), found no further case of long-term pe- 
riod change. Wood ( 1976) studied 45 stars and found small period 
changes with greater th an 99 per cent significa nce in 6 cases, but 
no large-scale evolution. iMattei & Fosteil 12000) performed a trend 
analysis of 383 long period variables using 90 years of AAVSO 
data. Only 9% showed evidence for trends, and of these only 9 
stars showed strong trends. IFercv & A3tl999) analysed a sample 
of 391 Mira variables and found some evidence for a slow average 
increase in periods with time, but few or no stars with large period 
variability. 

This suggests that over the past 50-100 years, period evolution 
has only occurred in ~ 1% of well-observed Miras. To find a case 
among the seven SC/CS stars classified in the GCVS as Miras with 
established periods is not ruled out by this rate. But the discovery of 
a second case implies an overrepresentation of the SC/CS stars. The 
rate of occurence of period evolution among semi-regular variables 
is not known. 

Period changes in Miras have norm ally been a ttributed to 
a recent thermal pulse or helium flash ( Woo d~& Zarrol 1798 lh . 
This explains the rarity of such objects. Fast period evolution 
will only occur during an d immediately after the helium flash 
IVassiliadis & Woodlll993l) . At any one time, only ~ 0.1% may 
be found at the rapid period increase. This model would therefore 
suggest that BH Cru is in a unique evolutionary phase, following a 
thermal pulse occuring during the 1960s. 

But the fact that, among the small number of SC stars, LX 
Cyg is also undergoing similar changes, makes the thermal-pulse 
explanation rather unlikely. The fact that historical records suggest 
recurrent period changes in BH Cru also argues against a thermal 
pulse. There are alternative models for period changes, predicting 
a non-linear instability in the Mira structure. In the case of R Hya 
I Ziilstra et al. 2002) those may yield a better fit to the data. 

In AGB stars, the visual pulsation amplitude is largely de- 
termined by varying molecular opacities: e.g., the extreme am- 
plitudes of oxygen-rich Miras is caused by the varying TiO and 
VO bands. The low molecular opacities in SC stars can explain 
their low pulsation amplitudes (Appendix|XJ. However, the molec- 
ular abundances become very sensitive to small changes in the 
chemical equilibrium, and this may affect the pulsation ampli- 
tude. An increased pulsation will lower the effective temperature 
and this shifts the chemical equilibrium. This process can result 
in a positive feedback, wh ich lea ds to larger pulsation amplitudes. 
Beddi ng! Conn & Ziilstral {2000) have suggested that a change in 
the amplitude can cause a period change (rather than the usually 
assumed reverse relation) due to non-linear effects. The positive 
feedback on the amplitude may thus affect the period as well, and 
may be a possible cause of the pulsation instability observed in BH 
Cru, LX Cyg and possibly VX Aql. 

Finally, although we argue against a thermal pulse, we note 
even if a thermal pulse did occur, the C/O ratio in BH Cru would not 
have changed. Carbon dredge-up following a thermal p ulse does 
not c ommence until ~ 250 yr after the thermal pulse iMowlavi 
1999), when the radius has returned to its pre-thermal pulse value, 
and is continuin g to decrease. The dredge-up for the models of 
iMowlavil \ 19991) lasts for about 100 yr while the period of the star 
is decreasing, well after the initial sharp period increase. Thus, in 
either case, the change from SC to CS star would be due to a lower 
temperature, and not a increase in the C/O ratio. The distinction be- 
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Figure 15. The wavelet analysis for LX Cyg. The lower three panels are on the same scale as Fig. [3 to facilitate comparison with BH Cm 



tween SC and CS stars may measure effe ctive temp erature, rather 
than an evolutionary progression l Stephenson 1973). 



6 CONCLUSIONS 

Our investigation of the evolving variable BH Cru shows an in- 
crease in its period of 25% within 25 yr. The period has stabilized 
at about 540 days. The visual semi-amplitude has increased simul- 
taneous with the period increase, up to a value of 1.25 mag. The 
changing period shows that the radius of the star has increased, and 
the temperature decreased, the latter confirmed by a slow redden- 
ing. 

The spectral type of BH Cru changed simultaneously from SC 
to CS. Chemical equilibrium modelling, both for the potosphere 
and a hydrostatic atmosphere, explains this as due to the decrease 
in the effective temperature. The lower temperature favours forma- 
tion of C2 and causes the fractional abundances of ZrO and YO to 
decrease. Our calculations explore a range of C/0> 1, for which 



oxides are still found although at reduced abundances. The distinc- 
tion between SC and CS stars does not require an evolution from 
C/0< 1 — >> 1, as sometimes is suggested, although this can 
also play a role. Infrared spectra suggest the possible presence of 
the iron-sulfide troilite in SC stars. Chemical calculations show that 
FeS is abundant in the photosphere, and in the absence of silicates 
or carbon grains, may be an important dust component. 

We determined new periods for a number of SC/CS stars, in- 
cluding three for stars with no previously known period. One star, 
LX Cyg, shows a much longer period than previously determined, 
and is confirmed as undergoing period evolution similar to BH Cru. 
VX Aql also shows an inconsistency between the present period 
and one previously reported. Among the few Miras known with 
evolving periods, BH Cru was unique in showing an increasing pe- 
riod. Its time scale for the evolution is ten times faster than that of 
the well-studied case of R Hya. With LX Cyg, a second case of 
rapid period increase is now known. 

Period changes in Miras are commonly attributed to thermal 
pulses. This appears unlikely in BH Cru, because other SC stars 
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show the same type of changes. As an alternative model we sug- 
gest the possiblity that a feedback between molecular opacities, 
pulsation amplitude and periods cause unstable periods among the 
long-period SC stars. 
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Table Al. GCVS data for the SC and CS stars . Sta rs in the last section 
are cla ssified as SC bv lCatchpole & Feasi 1 1 97 it) and lKeenan & BoeshaaJ 
119801) : the GCVS gives these as C stars. 
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a We note that the GCVS gives the wrong declination for this star. The 
correct values are a, S = 08 h 46 m 17.6 s , -71°02'20" (J2000) 
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Figure Al. Power spectrum (upper panel) of the CCD light curve (lower 
panel) of LQ Ara; the solid curve shows the best-fitting sinusoid. 
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APPENDIX A: PERIODS OF THE OTHER SC STARS 

The General Catalogue of Variable Stars (GCVS) lists only 1 1 SC 
stars (including BH C m) and 5 C S sta rs. These are listed in Ta- 
ble lATllKeenan & BoeshaaJ Jl98fj) and lCatchpole & Feastl <197ll) 
list a few other SC stars that the GCVS gives as carbon stars, and 
which perhaps should be classified as CS. In Table lAll these stars 
are in the third group (in italics). 

Of the 16 known SC and CS stars, 6 have no catalogued pe- 
riod. We attempted to determine some of their periods, and to con- 
firm existing periods, using new data (see below) and data from 
amateur databases. New light curves and/or power spectra were de- 
termined for the six stars, as follows. 



Figure A2. Power spectrum (upper panel) of the CCD light curve (lower 
panel) of AM Cen; the solid curve shows the best-fitting sinusoid. 

LQ Ara, AM Cen, UX Vol: Photometric observations of these 
three stars were obtained by one of us (VT) using CCD imaging 
with an unfiltered 100mm telephoto lens. For LQ Ara, the light 
curve and power spectrum (Fig. lAH indicate a period of 193 days. 
The amplitude is low, indicating that this star is a semiregular (not 
a Mira, as sua gested in the GCVS). For AM Cen (Fig.lA2l we find 
a possible period of 276 days, but this requires confirmation with 
more data. The star UX Vol (Fig. I A3 \ appears to be a semiregular 
with a period of 182 days. 

GP Ori: This star is classified as a semiregular (SRb) with a pe- 
riod of 370 d in the GCVS. We have analysed about 440 visual 
estimates from the AAVSO and VSOLJ databases, spread rather 
sparsely over the past 30 years. The power spectrum, shown in 
Fig. lA4l confirms GP Ori as a semiregular with two periods: 184 d 
(0.00543 d _1 ) and 339 d (0.00 295 d" 1 ). The pe riod ratio of 1.84 is 
typical of many semiregulars iKiss et allll99St) . The power spec- 
trum is somewhat complicated by yearly gaps in the data (the peak 
at 0.081 d _1 is an alias), and this is made worse by the fact that the 
shorter period is very close to half a year. N evertheless, both peri- 
ods seem secure. Bedding & Ziilstra ( 1998) showed that in double- 
period semiregulars, the longer period tends to agree with the Mira 
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Figure A3. Power spectrum (upper panel) of the CCD light curve (lower 
panel) of UX Vol; the solid curve shows the best-fitting sinusoid. 
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Figure A5. VSOLJ Light curve of VX Aql 

VX Aql is another example of an SC/CS star that has shown a large 
period increase. 

RZ Peg, S Lyr, R CMi: For each of these stars, visual obser- 
vations from AAVSO, VSOLJ and AFOEV databases show them 
to be Miras and confi rm the periods listed in the GCVS and in 
lMattei&Foste3 )2000). The data go back 80-90 years, over which 
time the periods are all stable. 
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Figure A4. Power spectrum of the light curve of GP Ori showing two 
distinct periods. The inset shows the spectral window. 



period-luminosity relation while the shorter period falls on th e 
semi-regular sequence (sequences C and B of IWood et alJll99^) . 
We also note that the power spectrum shows broadened envelopes 
that have been interpreted bv lBeddindiiool) as evidence for solar- 
like oscillations (i.e., stochastic excitation, presumably from con- 
vection). 

CY Cyg: We have analyzed 424 AAVSO observations from JD 
2438540 to JD 2452653. The star seems completely irregular, and 
the Lb classification (which is usually based on lack of data) seems 
warranted. 

VX Aql: This star is of particul ar interest because o f a C/O ratio 
suggested to be exactly unity lGreene&Wini3ll975l) . The GCVS 
lists it as a Mira but does not give a period. A VSOLJ light curve 
is shown in Fig. lA5l The AAVSO light curve can be found on the 
AAVSO web site. We have analyzed 86 observations from a sin- 
gle AAVSO observer which show the star is a Mira, with a pe- 
riod of 604 days and a very pronoun ced double maximum. Inter- 
estingly, Wvckoff & Wehin gei] (197^) adop ted a period of -~470 d 
which they attributed to Kurochkin 1 1958). It seems possible that 



